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Аппроксимация автоволновых решений в моделях ламинарного пламени” 
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Предложен новый алг оритм рас чёта параметров ламинарного пламени в г. орючей Смеси. От известных опти- 
мизационных алгоритмов разработку отличает способ аппроксимации точного решения задачи — постулиру- 
ется реалистичный вид не температурного поля, а его производной по автомодельной переменной. В ре- 
зультате процедура определения эффек тивных параметров модели упростилась, а погрешность аппроксима- 
ции сохранилась на приемлемом уровне. 
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Введение. Разработка действенных методов исследования автосолитонных уравнений представ- 
ляет большой интерес для ряда естественнонаучных и технических приложений [1]. Простейшей 
формой горения является плоское ламинарное пламя, формируемое процессами различной при- 
роды — химическими реакциями и микропереносом. Результирующий процесс представляет собой 
бегущую волну превращений [2]. Математической моделью пламени служат нелинейные уравне- 
ния с частными производными. В данной работе развит новый подход к решению уравнения го- 
рения, опирающийся на рациональный выбор функционального вида решения и отработанную 
технику оптимизации. 
Постановка задачи и описание моделей. Динамика автоволн горения описывается уравнени- 
ем [2]: 
2. 
РОТ +97 + шегет (1-т)’=0, Т(-®)=ь Г(®)=0, (1) 
02 07 

в котором температура Ги энергия активации ® выражены в единицах Т.„ (температуры сгорания 
в адиабатических условиях); л ^ 1 — эффективный порядок брутто-реакции; 2, И/и у — коэффи- 
циент диффузии, предельная скорость химического превращения и скорость распространения 
пламени соответственно. Решение (1) позволяет определять брутто-характеристики химизма пу- 
тём измерения скорости распространения пламени, его толщины и температуры. Определённые 
таким образом физико-химические характеристики ламинарного пламени используются при про- 
ектировании энергосиловых установок. 


Приблизим реальный химический источник, характеризующийся скоростью ИЁ“9” (1-7) =0, 
функцией ИТ? (1-Т): для удобства перейдём к безразмерным переменным Т =Т/Т.„ и 
е=е©/Т, 
смеси — своего рода форм-фактором источника, определяющим ширину благоприятного для ре- 
акции температурного интервала Тем + Тиах. Независимой переменной выбрана безразмерная 


температура Т, которая совпадает с безразмерной концентрацией [В] продукта В одностадийной 
реакции А > В+ О (О — тепловыделение). 


Последняя величина служит важной физико-химической характеристикой горючей 


ах " 
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Из условия совпадения максимумов исходной и приближающей функций следует связь 
модельного показателя В с содержательной величиной ©: 


29. 
в- 9/2. 1+(1+4/5) | (2) 
В отсутствие активации В —> 0 , при незначительной (по отношению к условиям в пламени) высо- 
те энергетического барьера реакции В -› ©"? , а в условиях, когда активация существенна, 


В-› ©. Безразмерная температура Т”, соответствующая максимальной скорости реакции в то- 
чечной системе — реакторе идеального смешения равна ВВ + 1) . Множитель а в приближающей 


_—е/т 


формуле определяется равенством аррениусовского е и степенного аТ? скоростных множи- 


телей при Т =Т”, из которого следует явное выражение: 


а= (1+1/Ве®. (3) 


0-1 


При высоком активационном барьере выражение (3) упрощается: а = е`?*". Однако практически 


интересен случай, когда величина © оказывается в пределах от одной до нескольких единиц. 
Модельный пример [3] с В=2 попадает в категорию реалистичных: ему соответствуют © =4 /3 


иа= (3 / 2е)' : 

Полученный таким способом коэффициент а следует включить в вероятностный фактор 
реакции И/ определяющий решение задачи (1). Поправка к скорости распространения пламени 
имеет химическую природу и оказывается существенной, если рассматривать горение моторных 
топлив. Например, для водородно-воздушных смесей величина ® изменяется от ^ 2 при стехио- 
метрическом соотношении компонентов до ^^ 15 на нижнем концентрационном пределе, что тре- 


бует корректировать скорость пламени. Соответствующий поправочный коэффициент близок к 
асимптотическому а? = е(®*? . 

Идея предлагаемого метода решения задачи состоит в том, что сначала выбирается фор- 
ма производной решения (1) по пространству, затем конструируется соответствующая функция 
источника, содержащая некоторые подгоночные параметры. Задача сводится к минимизации не- 
которого функционала совпадения (ФС) с применением стандартных численных методов. 

В качестве модельных решений — профилей установившейся температурной и (или) 
концентрационной волны — будем выделять: 1) квазилогистические, являющиеся решением 


уравнения аТ/а2 = сТ* (4-7), 2) тригонометрические, отвечающие уравнению аТ/а? = 
= са/псоз? [п(Г-1/2)] ‚ И 3) квазикинетические, соответствующие профилю температуры 


ат /а = се`®'Т (1-Г), где ©’ 0 и, как правило, меньше по абсолютной величине, у>1 — по- 
рядок реакции. 
Существенно, что во всех базовых моделях контур величины @Т/47 представляет собой 


пик, расположенный в интервале дозволенных значений Т ; на краях этого интервала реализу- 
ются однородные условия (начальные или конечные). Без потери общности можно считать, что 
профиль температуры и (или) концентрации исходных продуктов в волне горения задаётся фор- 


мулой вида сТ (Г)-(1-Т)', где Т (Г) — некоторая возрастающая функция, причём Т (0) =0. В 


последнем выражении учтено влияние энергетических и материальных факторов на структуру 
зоны активного горения. Выбранная математическая структура позволяет объединить все пере- 
численные профили волны горения, т. е. унифицировать рассмотрение. Явное выделение множи- 
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теля (1 № т) вызвано тем обстоятельством, что горение связано с материальными ограничениями 


и подчиняется реакции некоторого (не обязательно целого) порядка по реагентам. 
Рассмотрим волну горения с квазикинетическим профилем 


ат /а2 = ве-9" (1- т) (4) 


как вполне реалистичную. Вторая производная искомого распределения Т(2) по независимой пе- 
ременной равна 


ата? =в [6/7 -у/1-Т]е°" (и-Г)'|. (5) 


Подстановкой (4) и (5) в уравнение (1) нетрудно убедиться, что разрешающий его в выбранном 
классе профилей Т2) источник имеет вид 


геи [кс [в -+/(-7е°"(-7)] | (6) 


Зависимость (6) нуждается в определении подгоночных коэффициентов (, С, ©’ и у’. Их нетруд- 
но вычислить, зная вид источника: 


".„(Г)=И | ((- Г) (7) 

Минимизи ровав функционал совпадения 
ес - [["(Г)-=..(ГТот, (8) 

0 


с учётом необратимости реакции в пламени, которая определяется ограничением на параметры К 
и С обеспечивающим положительность множителя в фигурных скобках (6) при Т = 0 - 1, полу- 
чаем искомые параметры модельного источника {А, С, ©’, у’ }. Через них выражается скорость 
распространения пламени и — собственное значение краевой задачи (1): 

и=(ир/с)” к. (9) 
Полученный результат не связан с выбором профиля волны, а является следствием фундамен- 


тальных свойств уравнения (1) и регулярности предложенного способа его решения. Поскольку И/ 
и Р — заданные константы, определение скорости горения свелось нахождению нетривиального 


(К=1, С=0, 9'’„ 0, у’+у) аппроксимирующего выражения (6). 
Результаты компьютерного эксперимента. Численные эксперименты свидетельствуют о сла- 
бой зависимости эффективного порядка реакции в пламени от высоты активационного барьера 


О, что позволяет с хорошей точностью считать у’=у. Кроме того, оказывается, что ©’ близко 


к ©. Эти обстоятельства, с одной стороны, существенно облегчают определение подгоночных 
коэффициентов, а с другой — предоставляют основания для физической интерпретации расчёт- 
ных данных. В соответствии с (6) получается приближённое решение краевой задачи (1): 


ев" (1-РУ’ЧР - е+атяе, 
в= (ИС/р)”, 
0’ = 0,9390 - 0,0216?, т 
К = (1,028 — 0,0519) е-207146\°, 


с= е-!'23+6.0696+3,555у'" 
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Представление о погрешности решения (10) можно получить, сравнивая модельную функ- 
цию (6) при подстановке в неё коэффициентов (10) с исходной формой источника (7), а также 
само решение краевой задачи — скорость распространения пламени — с соответствующими ре- 
зультатами, полученными независимыми методами. Локальное отклонение предложенной мо- 
дельной функции от оригинальной составляет доли процента даже при наименее благоприятном 


сочетании параметров ©’ и у, когда максимальная скорость реакции реализуется при близкой к 0 


или, наоборот, к 1 безразмерной температуре 7’. В первом случае реализуются холодные пламе- 
на, во втором — нормальные или горячие. 

Сопоставлять непосредственно приближённое решение краевой задачи (10) с найденными 
альтернативными способами нецелесообразно, поскольку детали профиля источника (7), следо- 
вательно, качество их аппроксимации уравнением (6) и соотношениями (10) отражается, главным 
образом, на определяемом формулой (9) собственном значении и. Само же решение (профиль 
термохимических параметров в волне горения) является интегральной характеристикой, что сни- 
жает погрешность его определения. Поэтому представительно сравнивать вычисляемое нашим 
методом собственное значения и задачи с полученными известными способами. В безразмерном 
виде 

п =К/С!?. (11) 
Эта величина определялась многими авторами с использованием разнообразных методик [3]. 
Рассчитанная нами зависимость и =и (©, у) приближается формулой 


и (9 у) = е*778-18666* 2,032" а. 
: = . 
Удобными оказываются и менее точные приближения и = в (6, у): 
и (© у) — е?262-0,80596-0,446у ыы 
,у]|= 


ий (6, у) =: (1, 028— 0, а , (14) 


Формулы (12)—(14) хорошо согласуются с классическими результатами Зельдовича — Франк- 
Каменецкого, Кармана и другими приближениями [2—4]. 
Выводы. Результаты проведённых компьютерных экспериментов свидетельствуют о том, что по- 
лученные результаты приемлемо согласуются с классическими в физике горения. Расчётным пу- 
тём установлено, что сходное поведение демонстрируют все классы аппроксимирующих моде- 
лей — квазилогистические, тригонометрические и гибридные. Разнообразные модели и способы 
решения рассматриваемой краевой задачи приводят к близким скоростям ламинарного пламени. 
Этот факт свидетельствует о структурной устойчивости уравнения (1) и оправдывает применение 
предложенного нами нового способа решения задачи (1). 

Также установлены и (или) подтверждены следующие свойства решения задачи (1), (7): 

1) существует стационарный режим горения, которому отвечает некоторая скорость рас- 
пространения пламени и, полностью определяемая физико-химическими свойствами среды; 

2) среда неустойчива к параметрическим возмущениям, т. е. локальное возмущение уси- 
ливается и развивается в волну горения; 

3) пространственный масштаб волны нормального горения совпадает с диффузионной 
длиной, отвечающей времени химического превращения. 
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